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Abstract 

We use the theory of YFS resummation to compute the size of the expected resummed 
soft radiative threshold effects in precision studies of heavy particle production at the 
LHC, where accuracies of 1% are desired in some processes. We find that the soft QED 
threshold effects are at the level of 0.3% whereas the soft QCD threshold effects enter 
at the level of 20% and hence both must be controlled to be on the conservative side to 
achieve such goals. 
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The physics objectives of the LHC entail a precise knowledge of the Standard Model 
processes that are either background to the would-be discovery processes, such as Higgs 
and susy processes, or important for the normalization of the observed signal and back- 
grounds processes so that the cross sections of both types of processes can be determined 
in a way which maximizes the physics output of this pioneering high energy colliding beam 
device. One of the many sources of uncertainty that must be controlled are those then as- 
sociated with higher order radiative corrections to all aspects of the 15 TeV pp collisions. 
The dominant source of these effects are the higher order QCD corrections and many 
authors [1] have worked and are still working on these corrections. But at the level of pre- 
cision required for some of the LHC objectives, one must also check that the higher order 
EW corrections are under the appropriate control. It is well-known that initial state QED 
corrections can be large, particularly for resonate heavy particle production, and hence, 
already in Ref. [2,3] estimates were made of the size of the QED corrections to the struc- 
ture function evolution for the LHC energies. An effect at the level of a few per mille was 
found in most of the relevant range of the Bjorken variable x. More recently, in Ref. [4], 
a similar result was found with however a qualitatively different character - whereas the 
results in Refs. [2,3] all show that the QED correction decreases the structure functions in 
the relevant regime of x, those in Refs. [4] actually show a QED correction which changes 
sign in the relevant regime of x, but which is still similar in magnitude, a few per mille 
in general in the relevant regime. This could just be the difference in the bases used to 
present the respective results, as a recent more complete treatment of the these same 
effects in Ref. [5] confirms both of the behaviors found in Refs. [3,4]. In the current dis- 
cussion, we do not address these effects further. Rather, we focus on the potentially larger 
issue of the size of the non-universal QED corrections associated with the threshold be- 
havior of the heavy particle production at the LHC. To illustrate these effects, we treat as 
prototypical processes the processes pp — > V + n{ r f) +m(g) +X — > ££' + n'( r y) + m(g) + X, 
where V = W ± , Z,and £ = e,/i, I' — z/ e , v M (e, //) respectively for V = W + (Z), and 
£ = Vejffj,, £' — e,/i respectively for V = W~. These processes are potential luminosity 
processes for the LHC [6] 1 . What we want to do is to estimate the expected size of the 
threshold contributions from these corrections in the realistic acceptance for the process 
in the LHC detectors. In this way, we show what role if any these threshold radiative 
effects may play in precision LHC studies where the luminosity is expected to be needed 
at the few percent precision level. We recall for reference that the current precision on the 
FNAL luminosity is ~ 6% [7]. We stress that, for an overall precision of 2-3%, the theory 
budget in the error must be held to ~ 1% so that it does not figure too prominantly in 
the total LHC luminosity budget. 

Physically, the source of the effects which we calculate will be the following. In 
the complete calculation, assuming we have a factorization scale [ip comparable to the 
renormalization scale [ir, both of which are set near the heavy particle production scale 
M, one has the big logs L s = lnM 2 //^, L t h = L — lns/M 2 , where here we assume with 

1 We understand [7] that FNAL experiments are also considering these processes as a potential lumi- 
nosity processes. 
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Refs. [2,3] that the structure functions have been evolved from input data at scale yt, s that 
have QED radiative effects in them, s is the invariant mass squared of the hard parton 
production process. 2 Hence, all of the familiar big logs of the type In Q 2 /mj that one 
would expect from analogy with the big logs in e + e~ annihilation at high energies are 
removed by the factorization of the cross section's mass singularities as usual. Here, Q 
is any hard scale and m/ would be the respective quark mass. What do remain are the 
residual big logs L s , L t h, where the former are presumably all summed up by the structure 
function evolution. This leaves the latter, which we now study using the YFS methods 
in Ref. [9,10]. 

We start from the basic formula 

da exp (pp -> V + X — *■ It' + X') = / dxidxjF^XijF^x^da^XiXjs), (1) 

where the squared pp cms energy is s and the structure functions have their usual meaning 
here so that & exp (xiXjs) is the reduced YFS exponentiated V production cross section, V = 
W, Z. The YFS formula which we use in the latter cross section was already developed in 
Ref. [11], where the dominant higher order corrections from the QCD corrections for top 
production at FNAL were estimated. Here, we can extend that result by including the 
dominant QED corrections in addition to the ISR QCD corrections treated in Ref. [11]. 

We point out that the threshold resummation result which we obtained in Ref. [11] 
is agreement with the results of Catani et al. in Ref. [12] for the size of the higher order 
threshold resummation effects in QCD for tt production at the Tevatron. As we have 
explained in Ref. [13], this is not an accident. The Sterman-Catani-Trentadue (SCT) [14, 
15] resummation theory for QCD used in Ref. [12] is fully consistent with the QCD 
extension of the YFS theory that we discuss in this paper. As we explain in Ref. [13], 
one has to use the corresponding values of the j3 n residuals in the extension of the YFS 
theory to QCD and the corresponding higher order corrections to the infrared exponent 
SUM ir (QCD) in eqs. (15-18) in Ref. [13] to get an exact correspondence between the two 
resummation approaches. As we will see below, the result which we obtain for the QCD 
resummed threshold effects in Z production using our extension of the YFS methods to 
QCD will account for most of the actual QCD correction as given by the exact result 
through 0(a 2 s ) in Refs. [16,17]. This again shows the applicability of our methods to this 
problem and, by our correspondence in Ref. [13], the applicability of the SCT method as 
well. We do not pursue the latter method further here. 

Specifically, we recall the need to use Dokshitzer-Gribov-Lipatov-Altarelli-Parisi ( 
DGLAP ) synthesized [13] YFS formulas so that we avoid double counting of effects. In 
practice, this amounts to applying the factorization theorem to the YFS real and virtual 

2 We stress that there are apparently data available [8] from which the QED corrections from quarks 
have been removed so that the attendant QED correction to the structure function evolution would have 
a different boundary condition than that used in Refs. [2,3]. We are not aware of any results in the 
published literature where this type of calculation has been done. 
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infrared functions, removing those parts of these functions which generate the big logs 
L s . The factorization scheme used should be that used in the isolation of these big logs 
to the respective structure functions. This means that the functions B nls , B nls , S nls 
defined in Ref. [13] should be used in the YFS formulas applied to the factorized reduced 
hard parton-parton production cross section. Here, we stress that we are doing this for 
both the QCD and the EW QED corrections. This means that, in the basic YFS algebra 
illustrated in Ref. [11], we need to make the replacements, using an obvious notation, 

Tjnls r>nls , r>nls _ rynls 

D QCD ~ * D QCD D QED — QCEDi 

Tjnls rynls . r>nls rynls 

D QCD ~^ D QCD t n QED — D QCEDi 

onls onls I Qiils onls (c\\ 

D QCD ~^ '-'QCD ' '-'QED — '-'QCED \ z ) 

in the YFS exponentiation algebra in Ref. [11], with the attendant replacement of the 
reduced parton cross sections and structure functions to reflect the incoming pp state as 
it is compared to the incoming pp state in Ref. [11] and to reflect the V production and 
decay here versus the tt production there. 



In this way, we start from the basic result 

nr-. ™ n J «" 1 



n,m=0 ° j\ =1 



n 



.12-- 



f _^y_ e iy(p l + gi ~p 2 ~ q2 ^k n ~T,k , ]2 )+D QCED 



d 3 p 2 d 3 q 2 



i k' n J (2vr 



/3 n ,m{kl, • • • , k-n, ^l) • • • i k m ) n q , (3) 

P2 12 

where the new YFS residuals j3 n>m (ki, . . . , k n ; k[, . . . , k' m ), with n hard gluons and m hard 
photons, represent the successive application of the YFS expansion first for QCD as in 
Ref. [11, 18, 19] followed by that for QED as some of us have used in many applications 
as in Refs. [10]. The infrared functions are now given by 

SUM IR (QCED) = 2a.toBgZ ED + 2a a B$ ED 

/dk 
- (e-^ - 9(K max - k )) S$. ED (4) 

where the dummy parameter K max is just that - nothing depends on it. We have taken 
the liberty, not to be viewed as permanent, to set the dummies in QCD and QED to be 
the same value. The infrared algebra realized here for QCD and QED we will sometimes 
denote as Quantum ChromoElectroDynamics (QCED). 

The result J3J) may be realized in the context of (0) using Monte Carlo methods as 
illustrated in Refs. [19]. We shall present those methods in detail elsewhere. Here, to 
illustrate the size of the QED effects in the presence of the QCD resummation as well, we 
use a semi-analytical evaluation of ((TJ. 
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Specifically, in order isolate explicitly ISR effects, we take the case of single Z produc- 
tion and refer the case of single W production to Ref. [20] . We also call attention to the 
complete EW 0(a) corrections to single heavy gauge boson production in hadron collid- 
ers in Refs. [21-25]. We work in the approximation that, when we retain the dominant 

infrared effects at the /?o°o° -level, these two sets of corrections commute with one another: 
physically, gluons are EW singlets and photons are color singlets. Thus, the only issue is 
whether particles such as quarks with both color and EW charge have significant higher 
order corrections in which the two sets of corrections do not commute. Here, we appeal to 
the Heisenberg uncertainty principle and note that the average soft photon energy fraction 
in the YFS formalism is x avg (QED) = 1 (QED)/{l+^{QED)) and the average YFS soft 
gluon energy fraction is correspondingly x avg (QCD) = j(QCD)/ (l+ r y(QCD)), where the 
YFS radiation probability strengths are given by 7(A) = ^^(L s - 1), A = QED, QCD. 
For quarks, Ca = Q}, Cf, respectively, for A = QED, QCD, where Qf is the quark elec- 
tric charge and Cf is the quadratic Casimir invariant for the quark color representation. 
Here, we variously use (Xqed = ca, &qcd = Oi s . As the ratio of the two fine structure 
constants is ~ 10 in the LHC/Tevatron environment in which we work, we expect that 
the typical time for the dominant QCD corrections we exponentiate to occur is more 
than an order of magnitude shorter than the corresponding time for the analogous QED 
corrections which we exponentiate. Thus, we do not expect that significant effects in our 
leading exponentiated analysis are missing due to sets of QCD and EW QED corrections 
that do not commute. We then compute the ratio r exp = a exp / ' o~Bom with and without the 
QED contributions, where, for definiteness, we note that the partonic Born cross section 
we use here is well-known and we generate it by the standard methods from eq.(43) in 
Ref. [26], with the understanding that we retain only the Z contribution and that we 
substitute the respective incoming quark chiral couplings for the incoming e— chiral cou- 
plings accordingly. (We stress that we do not use the narrow resonance approximation 
here.) We get the results 

'1.1901 , QCED = QCD+QED, LHC 
1.1872 ,QCD, LHC 

1.1911 , QCED = QCD+QED, Tevatron ^ 
1.1879 ,QCD, Tevatron 

We see that indeed the QED threshold correction is at the level of 0.3% at both the 
Tevatron and at the LHC and is in principle significant for studies at the percent level such 
as those that are envisioned for the precision luminosity determination for the LHC [6] 3 . 
We also corroborate here the well-known [16,17] result that the QCD threshold correction 
is needed through 0(a 2 s ) for such studies. The results in © are entirely consistent with 
the exact 0(a 2 s ) results in Ref. [16,17], where comparison with the pure QCD results in 
eq.(0) shows that our threshold effects reproduce most of the complete exact results at 

3 We checked that the variation of the renormalization scale, which we set equal to the factorization 
scale, above and below its nominal value Mz by a few widths of the Z does not change our result that 
the QED ISR threshold correction enters at the .3% level 



f exp 



4 



both the LHC and the Tevatron. The size of our exponentiated QED effects in (JHJ) is fully 
consistent with the exact O(a) results in Refs. [21-23], where the ISR QED correction was 
found to be 0.4%. In arriving at our results in JHJ), we have used the parton distributions 
from Ref. [27]. 

Our result for the size of the QED threshold effect is similar to the size found in 
Refs. [2-5] for the size of the QED effects on the structure function evolution itself. From 
these two results, we can conclude that the ( higher order ) QED correction to heavy 
particle production at the LHC will be at the several per mille level and it can not really 
be ignored for percent level precision studies. 

In conclusion, we have introduced a new theoretical construct in this paper, QCED, 
the regime of the large IR effects in the large momentum transfer regime, where we expect 
the QED and QCD exponentiation algebras to commute to a large extent and to the 
extent that they do we identify the resulting theory as Quantum ChromoElectroDynamics. 
Corrections to it are then to be treated perturbatively in the respective combined YFS 
residuals Pnm ED - We have used this paradigm to estimate the size of the QED threshold 
effects in single gauge boson production at the LHC type energies. We find that such 
effects are large enough that they must be taken into account for precision studies, such 
as those envisioned for precision luminosity determinations. 
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